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Abstract 

Matter density uncertainties can affect the measurements of tlie neutrino oscillation pa- 
rameters at future neutrino factory experiments, such as the measurements of the mixing 
parameters 6'i3 and 6cp- We compare different matter density uncertainty models and dis- 
cuss the possibility to include the matter density uncertainties in a complete statistical 
analysis. Furthermore, we systematically study in which measurements and where in the 
parameter space matter density uncertainties are most relevant. We illustrate this discus- 
sion with examples that show the effects as functions of different magnitudes of the matter 
density uncertainties. We find that matter density uncertainties are especially relevant for 
large sin^2^i3 > 10~^. Within the KamLAND-allowed range, they are most relevant for the 
precision measurements of sin^ 2^13 and 5cp, but less relevant for "binary" measurements, 
such as for the sign of Am'^i, the sensitivity to sin^2^i3, or the sensitivity to maximal 
CP violation. In addition, we demonstrate that knowing the matter density along a specific 
baseline better than to about 1% precision means that all measurements will become almost 
independent of the matter density uncertainties. 



''Email: tommy@thcophys.kth.se 
''Email: wwinter@ph.tum.de 



1 Introduction 



Neutrino oscillation physics has entered the era of precision measurements of the leading 
atmospheric {Arri^^ and ^23) and solar {Am^i and 612) neutrino oscillation parameters. 
Especially, the recent results of the Super-Kamiokande [I1I21, SNO Oil], and KamLAND j3] 
experiments have introduced this era. In addition, the leptonic mixing angle ^13 has turned 
out to be small by the results of the CHOOZ experiment [HllZl- Currently operating or 
planned conventional beam experiments [SliniCni and planned superbeam experiments fH] 
IT^IT!^ will complement this information by further measurements of 6'i3, the neutrino mass 
hierarchy, and the leptonic CP violating phase 6cp- Finally, future neutrino factories (for a 
summary, see ^1]) could find ^^13 > and test leptonic CP violation even below sin^ 2^13 < 

10"^ [ig. 

Because of the high precision of these future experiments, it is important to study the im- 
pact of matter density effects on neutrino oscillations, since they may affect the different 
measurements in a substantial way. It has for a long time been known that neutrino oscil- 
lations are influenced by the presence of matter |16 | IT7 | ITH]. and therefore, matter density 
effects on neutrino oscillations in the Earth have been thoroughly investigated in different 
contexts and with a variety of models (see, for example, Ref. and references therein). 
In most of these analyses, the Preliminary Reference Earth Model (PREM) matter density 
profile [201 been used, which has been obtained from geophysical seismic wave mea- 
surements. However, the knowledge on the PREM matter density profile along a certain 
baseline through the Earth's mantle is limited from geophysics to about 5% precision PHl22j . 
which means that matter density uncertainties can affect the precision measurements of the 
neutrino oscillation parameters, such as the leptonic mixing parameters {612, On, ^23, and 
^cp) as well as the neutrino mass squared differences (Am^i and Am|]^). 

In this work, we will mainly focus on the effects of the matter density uncertainties on the 
future measurements of 613, sgn(Am|]^), and 6cp, since these measurements are the most 
interesting ones for future long baseline experiments. Especially, neutrino factories could be 
affected by matter densities for two reasons. First, they will be operating in the statistics 
dominated regime with a very high precision. Second, in comparison with superbeams, they 
are often proposed with very long baselines L > 700 km with energy spectra covering the 
matter density resonance in the Earth's mantle. Thus, they will be strongly influenced by 
matter density effects themselves (see, for example, Refs. pH|l24p25p2(jp27p28p29p30p3ip32p33j ) 
and also by matter density uncertainty effects |Mll^EglE7IIISEHilMlllIlliIlll21llgllIglllllll^ 

uni- 

The paper is organized as follows. In Sec. |21 we classify the in the literature existing models 
for matter density uncertainties into three different categories. Next, in Sec. El we give 
the experimental description of the neutrino factory setup that we have been using for our 
analysis as well as we describe how the simulations have been carried out. Then, in Sec. |31 
we present a qualitative discussion of the matter density uncertainty effects in the neutrino 
oscillation transition probabilities for the appearance channels z/g — > t'^, z/^ — > z/g, z/g — * z/^, 
and z/^ — > Ue- In Sec. we study the impact of matter density uncertainties on the most 
interesting future measurements at a neutrino factory. Furthermore, we show examples of 
such measurements. Finally, in Sec. (HI we present a summary of the obtained results as well 
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Figure 1: The model of the average matter density as a measured quantity, illustrated for 
a baseline of L = 7 500 km and a matter density uncertainty of Ap* = 5%- p. The deviation 
from the average matter density p, i.e., the quantity Ap, is measured by the experiment 
within the gray-shaded area. This area is determined by Ap* = p (1 ± 5%) at the la confi- 
dence level with Gaussian statistics. A possible profile from the random fluctuations model 
in Ref. HV^ is shown as the dashed curve. 



as our conclusions. 



2 Modeling matter density uncertainties 

In most experimental simulations, the matter density is assumed to be constant along the 
baseline: First, taking N matter density layers instead of one essentially increases the 
computational effort by a factor of N. Second, the differences in the results are, especially for 
not too long baselines, minor compared to the increased effort and the effects of correlations 
and degeneracies. Therefore, the matter density profile effect, i.e., the difference between 
a constant and a realistic matter density profile, is often not taken into account. However, 
for a fit of the data for a future long-baseline experiment it should not be a problem to 
incorporate it at least at the level of the PREM matter density profile. 

In this work, we are interested in the uncertainties of such a profile, which are coming from 
geophysics. These uncertainties have different effects than the matter density profile itself 
and they would essentially act as systematical errors in the measurement. However, the 
matter density profile can, up to a certain level, also be absorbed as such an uncertainty 
into the constant matter density |12]- In the following sections, we will therefore mainly 
study the consequences of matter density uncertainties as perturbance to a constant matter 
density profile, which should for specific experiments, of course, be replaced by the profile, 
which is as close to the realistic profile as possible. 

In the literature, there exist, in principle, three different classes of models for matter density 
uncertainties: 
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Statistical models to study the effects of realistic fluctuations. These models are av- 
eraging over the effects of many profiles with similar properties to be as realistic as 
possible. Then, the average characteristic effects of matter density uncertainties on the 
neutrino oscillation transition probabilities or other quantities of interest are studied. 
Examples are the path integral method to average over many profiles in Refs. P7|l39p44j 
and the random fluctuations model in Refs. fOllS]- In this work, we use the latter 
one, which demonstrates the dependence on the length scale A and the amplitude 
Ap of the fluctuations. In this approach, the actual value of the matter density p 
is fluctuating around the average matter density p with the length scale A and the 
amplitude Ap following some (truncated) Gaussian distributions with standard de- 
viations a\ and ciAp, as it is illustrated in Fig. [T] with the dashed curve. From the 
discrepancies of realistic geophysics proflle reconstructions 21J, one can estimate a set 
of standard parameters for these four parameters, which we will use if not otherwise 
stated: A ~ 2 000 km, ~ 3/4 A, Apo ~ 3%p, and a/^p ~ 1/3 Ap. 

Simulation models assuming the matter density to be unknown. The matter den- 
sity proflle is assumed to be unknown within certain limits in order to simulate the 
matter density uncertainty effects on neutrino oscillation parameter measurements. 
In these models, the unknown matter density proflle (up to a certain degree) im- 
plies new parameters to be determined by the measurements, which, of course, are as 
well correlated amongst themselves as they are correlated with the neutrino oscilla- 
tion parameters to be measured. Examples include the Fourier expansion method in 
Refs. |38|I3U1 E^ and models assuming an uncertainty for the average matter density 
in Refs. j^EOlllSllISlllS] . In Ref . jl2] , it is shown that the most interesting effects of 
matter density uncertainties can be translated into an uncertainty (or a shift) of the 
average matter density along the baseline. In fact, the matter density proflle effect 
can also be incorporated into the model by this approach, where an uncertainty of 5% 
on the average matter density is quite a safe choice for not too long baselines. The 
"model of the measured mean matter density" , which is illustrated in Fig. Q uses the 
average matter density from the PREM matter density proflle and allows it to vary 
within some range Ap*. It is chosen to be Ap* = 5% ■ p in Refs. I^OlllHllISlllSl using 
this model in order to safely cover the effects of the matter density uncertainties and 
even the proflle effect for not too long baselines. This means that the average matter 
density is measured by the experiment within 5% externally imposed precision. From 
Fig. ^ it should be clear that the random fluctuations can also be simulated with 
this model, since they are partially averaging out. However, it can be shown that for 
too large interference effects among the different matter density layers, such as for 
baselines L > 5 000 km, the shape of the spectrum changes too much to simply ignore 
the matter density proflle effect. Thus, experiments must, in their data analyses, take 
into account the proflle effects for such long baselines even within this model. 

Models for a specific baseline. For a speciflc baseline, one can try to obtain further 
geophysical information in order to incorporate the matter density proflle as good as 
possible (see, for example, Ref. [IE])- The deviations from the PREM matter density 
proflle can in this case be rather large. However, the matter density uncertainties of 
such a proflle are drastically reduced. 
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Figure 2: The mapping between the length scale X of the random fluctuations model 
and the quantity Ap*/p of the measured mean model for different fluctuation amplitudes 
and baselines. For the other parameters from the random fluctuations model in Ref. liy^ 
the standard parameters therein were used: ax = 3/4 A and a^p = 1/3 Ap. The mapping is 
performed by minimizing the -function for the NuFact-II experiment (cf., Sec.\^ between 
the energy spectrum of a random matter density profile and an average matter density profile 
with the density p + Ap* with respect to the matter density shift Ap* . It turns out that the 
minimal Ax^ is, for short enough baselines, much smaller than unity, which means that the 
random matter density profile can be approximated by an average matter density profile to 
a very good approximation. The shown functions are computed by averaging the shifts Ap* 
over 2 000 such random profiles and using interpolation. 

From the above discussion, it should be obvious that the models within each class are rather 
similar to each other, although they may be using different approaches and are aiming at 
different goals. However, one can also transfer knowledge from the statistical models in order 
to use it for a better estimation of the parameters in the simulation models with unknown 
matter density parameters. For example, the random fluctuations model (which simulates 
the effects of the geophysical matter density uncertainties for different length scales and 
amplitudes to be obtained from geophysics maps) can be used to estimate the precision of 
the matter density within which it is assumed to be unknown. Let us assume that we can 
identify the fluctuations to have a certain length scale A and amplitude Ap, what is the 
deviation of the average matter density Ap*, which needs to be used in the simulations? 

As it is noted in Ref. jl2], the energy spectrum with the fluctuations can, for short enough 
baselines, be fit to the one with a constant, but different, average matter density. The 
average difference between the original and fit average matter densities can then be taken 
as a measure for the allowed error on the average matter density Ap*. In Fig.|21 this mapping 
is performed as a function of the length scale A for two different baselines and two different 
amplitudes of the fluctuations. In this figure, it can clearly be seen that the fluctuations 
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partly average out for short length scales. However, for length scales X ^ L there is no 
difference between the random fluctuations and the measured average models, since they are 
identical in this limit. It is also demonstrated that in many cases a somewhat smaller value 
of Ap* than 5% should be sufficient for a fluctuation amplitude of Ap* = 5% ■ p. However, 
in this case, the matter density profile effect is not taken into account and it should be kept 
in mind that we are only comparing average statements. 

Closing this argumentation, the different length scales A also correspond to the length scales 
of certain Fourier modes, which means that fluctuations at a certain length scale leads to 
an uncertainty of certain Fourier coefficients. These can by their correlation (see Ref. |12]) 
be translated into an uncertainty of the leading Fourier coefficient. This means that the 
strength of this correlation between the higher order Fourier coefficients and the leading 
ones is another qualitative interpretation of Fig. |21 We will therefore further on focus on the 
effects of the average matter density as a measured quantity in a full statistical analysis, since 
this should simulate the matter density uncertainties very well and it is suitable for more 
sophisticated statistical simulations including correlations and degeneracies (c/.. Sec. EI). 

3 Description of experiment and simulation 

Although superbeams are often proposed with quite long baselines, they are either operating 
in the statistics dominated regime (first-generation experiments) or are limited by the back- 
grounds and their uncertainties (superbeam upgrades). Therefore, the class of experiments 
for which matter density uncertainties are most interesting is neutrino factories: First, they 
usually have long enough baselines, which means that matter effects substantially influence 
the appearance rates. Second, their precision is very good, which means that matter density 
uncertainties become a major impact factor. In order to demonstrate the effects of matter 
density uncertainties, we use the advanced stage neutrino factory setup NuFact-II from 
Ref. |in| with a muon energy of 50 GeV. This setup has a target power of 4 MW, correspond- 
ing to 5.3 ■ 10^° useful muon decays per year. In addition, we assume eight years of total 
running time, four of these with a neutrino beam and the other four with an antineutrino 
beam. For the detector, we simulate a magnetized iron detector with an overall fiducial 
mass of 50 kt. The baseline length is assumed to be L = 3 000 km, which has turned out 
to be reasonable for CP measurements. We use the analysis technique, which is described 
in the Appendices A, B, and C of Ref. ^D], including the beam and detector simulations. 
As described above, we use the average matter density p = 3.5g/cm^, which corresponds to 
this baseline, and we allow an uncertainty of Ap* = 5% ■ p. Later, we will discuss the effect 
of this uncertainty in greater detail. This means that the matter density is treated as all the 
other neutrino oscillation parameters to be measured by the neutrino factory. However, the 
average matter density is assumed to be known from external measurements within some 
uncertainty, and therefore, the matter density is considered as an external input coming 
from geophysics (see also Fig. Q). This treatment does not only take into account that the 
matter density is not precisely known, but it also allows correlations between the neutrino 
oscillation parameters and the matter density (c/., Ref. |12])- Since we include all relevant 
appearance and disappearance channels in our analysis, the only parameters which are mea- 
sured much better with a different class of experiments are the solar parameters. Hence, 
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we use a similar technique to include the external input from the solar parameters and we 
assume that their product, which the neutrino factories are sensitive to, is measured with 
15% precision by the KamLAND experiment by the time the neutrino factory experiment 
would take place |T71 H5]. 

For the analysis of the quantities of interest we include, where applicable, systematics, multi- 
parameter correlations, and degeneracies as it is described in Ref. [lU] for the individual mea- 
surements. To summarize, the treatment of systematics is done in a straightforward way by 
minimizing the x^-function over the auxiliary systematic parameters. Multi-parameter cor- 
relations are included by projecting the n-dimensional fit manifold onto the axis of interest. 
This approach always needs to be applied if one wants to know the precision of an individual 
parameter, since the other parameters (such as the CP phase) cannot be determined better 
by the experiment itself and all possible points within the fit manifold are equally allowed.^ 
Alternatively, one could, as a matter of taste, also project onto the plane of two parame- 
ters in order to have the information less condensed, but this approach makes it harder to 
discuss the dependence on the other parameters, such as the true parameter values. As an 
important difference to the analysis of existing experiments, the reference data set (event 
rate spectrum) of future experiments is generated with a given set of true parameter values 
assumed to be realized by Nature and it is not provided by a measurement. In many cases, 
the results strongly depend on the true parameter values, which can vary within their cur- 
rently allowed ranges. It is therefore our philosophy to condense the information as much 
as possible in order to be able to discuss the true parameter value dependencies, which are 
especially relevant to systematically assess the potential of future experiments. 

As far as degeneracies are concerned, we know the (^cp, 6*13) jSS], sgn(Am|]^) and 
(6*23, '?r/2 — 623) ISO] degeneracies, which lead to an overall "eight-fold" degeneracy 
Since the current best-fit value for the atmospheric mixing angle is 623 = vr/4, we only 
consider the first two degeneracies. The definition of the quantity of interest quite often 
already includes the way the degeneracies have to be treated. For example, defining the 
sensitivity to sin^ 2^13 as the largest value of sin^ 26*13, which cannot be distinguished from 
sin^ 2^13 = (at the chosen confidence level) already implies that any value of sin^ 2^13 
fitting sin^ 2^13 = of any degeneracy is below this sensitivity limit. Thus, for sin^ 26*13 
above this sensitivity limit it is guaranteed that we will find sin^ 26*13 > with the proposed 
experiment at the chosen confidence level. However, in some cases, including degeneracies 
is not that straightforward. For example, for the precision of the measurement of sin^ 2^13 
we do not find a simple way to include them, and therefore, we only show the results for 
the best-fit manifold. Similar to the correlations, an existing experiment would analyze its 
data in a different way. Such as it was performed for a long time with the solar parameters, 
different isolated islands in the parameter space would be allowed and subsequently reduced 
by the future experiments (or the combination of some experiments). However, it is the goal 
of the analyses of future experiments to keep the number (determined by the degeneracies) 
and the extensions (determined by the correlations) of these islands as small as possible by 
optimizing the experiments before they are built. Therefore, it is reasonable to condense the 

^If the other parameters (that the quantity of interest is correlated with) were better determined by 
preceding measurements, then they should be either included as external input or, even better, by using 
the x^-function of the relevant experiments. This does especially not apply to the CP phase Sqp- 
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information as much as possible in order to discuss the dependencies on the true parameter 
values. 

For the neutrino oscillation parameters, we choose (if not otherwise stated) the best-fit values 
Amj^ = 3.0 • 10-3 eV^ sin^ 2^23 = 1.0 [S21, Am|i = 7.0 ■ IQ-^eV^, sin^ 2^12 = 0.8 (see, for 
example, Ref. ^^), as well as we only allow values of sin^ 26'i3 below the CHOOZ bound 
and we assume a normal mass hierarchy, i.e., sgn(Am|^) = 1. Because of the symmetric 
operation of the neutrino factory, results would not differ much qualitatively for an inverted 
mass hierarchy. In general, we do not make any special assumptions about the true value 
of the CP phase Sep- However, in some figures, we will choose certain values for illustrative 
purposes, since studying the dependence on the CP phase is not the subject of this work. 



4 Qualitative discussion of the appearance channels 



The goal of this section is to arrive at a qualitative understanding of the effects of matter 
density uncertainties. However, we will see that there are quite many factors involved in 
this issue, which we will discuss in the following in somewhat greater detail. At the end of 
this section, we then will conclude with some very basic qualitative observations. 

For long-baseline experiments, the appearance probability Papp = Pefi, Pfj.e, Psp, or P^g in 
matter can be expanded in the small mass hierarchy parameter a = Am'^i/Am'^i and the 
small leptonic mixing angle 6^13 up to second order in these parameters as 

Papp = Pi±P2 + P3 + P4 

~ .1.2 9^ Sin^[(l--4)A] 

~ sm /t7i 3 sm (793 

(1 - Ay 

A sin(iA) sin[(l - i)A] 

± a sin 26*13 sin 26*12 cos 6*13 sin 26*23 sin ocp sin A 

A (l-A) 

. . n . r X Sin(iA) sin[(l - A) A] 

+ a sm 26^13 sm 26^12 cos 6^13 sm 26^23 cos Onp cos A 

A (l-A) 



+ c? sin^ 26'i2 cos^ 6'23 



sin' (A A) 
I2 ' 



Here A = Aml^L/{AE) and A = ±{2^G pUeE) / Aml-^ with Gp being the Fermi coupling 
constant, L the baseline length, E the neutrino energy, and Uf, the (constant) electron 
density in matter. This expansion is qualitatively valid independently of the channel used, 
but the signs before the second term and in the definition of the parameter A depend on the 
direction z/g t'^ and if one is using neutrinos or antineutrinos. The sign before the second 
term is positive for z/g and neutrinos or — ^ z/g and antineutrinos and it is negative for 
^'/i ^^'-^ antineutrinos or —>■ v^. and neutrinos. The sign in A is positive for neutrinos 
and negative for antineutrinos. Using parameter values for a within the KamLAND-allowed 
range, it turns out that for long-baseline experiments all four terms in Eq. contribute to 
the appearance channels with similar orders of magnitude. However, the relative weight of 
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these four terms is determined by the true values of the parameters a and ^13. 

In Eq. (PJ, the matter potential enters via A and the matter density resonance condition can 
be expressed as A ~ 1. At the resonance, the factor sin[(l — y4)A]/(l — A) becomes maximal, 
which means that the first three terms of Eq. (Q) are strongly influenced by the resonance 
condition. Especially, the quadratic dependency in the first term of Eq. ([T]) makes it peak 
much stronger at the resonance than the second and third terms, and therefore suggests 
that it is most sensitive to a change in any of the parameters in A. The effect of such a 
parameter change would also create a shift of the position of the resonance in the neutrino 
energy spectrum. This is, for example, often used for the sign of Amg^^-measurement, since 
A changes sign for the inverted mass hierarchy, producing a very distinctive solution. In 
this case, the resonance is essentially shifted towards negative neutrino energies, which is, 
of course, an unphysical solution. Since the first term of Eq. is proportional to sin^ 2^13, 
mass hierarchy measurements are therefore favored for large values of sin^ 26'i3 (and small 
values of a in order to keep the other terms small, which cause correlation and degeneracy 
problems for this measurement). Another parameter, which enters in A, is the electron 
density Uf. representing the matter density p? Similar to the sign of Amg^^, we therefore 
expect that the neutrino energy spectrum is especially deformed by the first term in Eq. (^J) 
for a shift of the matter density. 

In order to demonstrate the above, let us investigate the impact of a perturbation (or shift) 
IS.A on the parameter A somewhat closer. This shift corresponds to the amplitude of 
the matter density uncertainty Ap. One can write the absolute shift of the appearance 
probability Papp up to first order in the relative shift of the parameter A as:^ 

APapp = Papp(i + Ai)-Papp(i) 

2 ( AAcot[(l - i)A] - ^-J^ ^1 

1 — 1A\ 

+ ( -AA cot[(l - A)A] + AA cot(AA) - ^ — -J- \ P2,3 



+ 2 [aK cot(iA) - 1) P4 ^ + ((Ai/i)2) , (2) 



where P2^3 = ±P2 + P3. The absolute error of the appearance probability Papp is here an 
appropriate representation for the statistical analysis, since the change in the event rates 
and their impact on the x^-value is proportional to the absolute (and not relative) error. 
Note that all prefactors appearing together with the Pj's [i = 1, 2, 3, 4) in Eq. © are nearly 
of the same order of magnitude, at least for neutrino energies E > 20 GeV, as it can be 
seen from Fig. 121 However, the orders of magnitude of the different terms in Eq. (0) can 
still be very different from each other, since these terms are proportional to the absolute 

^The electron density rig is related to the matter density p as = -^f^P, where Y^, is the average number 
of electrons per nucleon (in the Earth: Ye ~ \) and rriM ~ 940 MeV is the nucleon mass. 

•^A comparison with an exact numerical calculation suggests that this first order approximation is sufh- 
cient for our purposes. 
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Figure 3: The pref actors of Pi, P2,3; and in Eq. |^ as functions of the neutrino energy 
E fSGeV < E < 50GeV; for the parameter values Amii = 3.0 ■ 10-3eV^ L = 3 000 km, 
and p = 3.5g/cm'^. The vertical gray line marks the value of the resonant neutrino energy 
Eres ^ 11 GeV. 

values of the Pj's defined in Eq. which per se could be very different. For example, 
for large values of sin^ 26'i3 the first term in Eq. is relatively large and would therefore 
also be strongly affected by matter density uncertainties. In addition, all terms in Eq. ^ 
are proportional to A A/ A to first order, which means that the absolute change of one of 
the terms is obtained by multiplying the prefactor with the relative change of the matter 
density and the corresponding Pj itself. In Fig. El we observe that all prefactors are finite 
and not singular at the resonant neutrino energy E^^s — 11 GeV and that the prefactor of 
P4 is a constant with respect to the neutrino energy E. Most interesting, the prefactor of 
Pi changes sign at the resonance, which means that a small shift AA of the matter density 
close to the resonance will cause 

• a spectral distortion and a change of the normalization in the first (and also somewhat 
in the second and third) term 

• only a change of the normalization in the fourth term 

of Eq. (P). Since the effects of the matter density shift can be of equal size, this observation 
makes the analysis of matter density uncertainties rather complicated. However, it should 
be noted that especially the first term can also simulate spectral distortion effects in addition 
to a change of the normalization, which means that it allows more sophisticated correlations 
than the other terms. 



9 



As far as correlations are concerned, we observe that in Eq. (P) only the combinations 

= s.n2..3 ^'°'"-f^' . (3) 

1- A 

^ . ^„ sin(AA) 

C2 = a sm 26^12 — — - (4) 

A 

are present. In Eq. (Q), Cf appears in the first term, C1C2 in the second and third terms, 
and C| in the fourth term. This means that a change in any parameter value in these 
terms can be compensated by another one. Since sin^ 2^13 is the quantity of interest to 
us, it will lead to correlations with the matter density shift AA. However, the quantities 
a and sin 2^12 are given by external measurements with lower precisions than the matter 
density, which means that they can easily absorb a change of the matter density from the 
correlation in the first term of Eq. (0). Therefore, especially for large sin^2^i3, where the 
influence of matter density uncertainties to the first term in Eq. (jT)) is large (c/., Eq. 0), 
the combination Ci can have a major impact on measurements. 

Another important ingredient in the analysis is statistics. This means that the impact of 
the matter density uncertainties can be suppressed in statistics dominated regions of the 
parameter space. Since some measurements are limited by statistics, the matter density 
uncertainty will not at all have a strong impact on these measurements. For example, the 
sin^ 2^13-sensitivity limit describes the ability for an experiment to establish sin^ 2^13 > 0. In 
this case, the reference rate vector is computed with sin^ 26*13 = and the statistics from the 
first three terms in Eq. which are only present for sin^ 26*13 > 0, limits the measurement. 
Thus, the sin^ 26'13-sensitivity limit is constrained by the absolute size of these terms, which 
cannot be compensated by the correlation with matter density uncertainties because of the 
spectral dependence of the first prefactor in Eq. (0). Since all of the quantities of interest to 
us are suppressed by sin 26*13 in Eq- ©, a similar argument can be used for small values of 
sin^ 2^13, where most measurements are limited by the statistics in the appearance channels 
and not by the matter density uncertainties. In addition, matter density uncertainties 
will turn out not to be relevant for many measurements with a strong "binary nature", 
i.e., measurements which only have two very distinctive different answers. The ability to 
distinguish these two answers often depends more on statistics than on the matter density 
uncertainty, which cannot simulate the different solution. One such example is the mass 
hierarchy determination discussed above, which could only be mimicked by matter density 
uncertainties larger than 100% ■ p. 

Putting all pieces together, we have found that especially sin^ 26*13 as a quantity of interest 
is highly correlated with the matter density uncertainties. However, the effects of the 
matter density uncertainties are proportional to the four terms in Eq. themselves, which 
means that they are (for the quantities of interest) suppressed by sin^ 26*13 itself. Therefore, 
we expect a large influence of matter density uncertainties for large sin^ 26*13 and a small 
influence for small sin^ 2^13 because of domination of statistics. For the sin^ 2^13-precision 
measurement the correlation between sin^ 2^^i3 and the matter density uncertainties should 
be limiting the measurements for large values of sin^2^i3. However, for CP measurements 
the first term in Eq. acts as a background. The matter density uncertainties should 
introduce an additional systematical error for large values of sin^ 26*13, such as a background 
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uncertainty. In addition, 6cp is highly correlated with sin^ 2^13, which again is correlated 
with AA, which suggests that it is also indirectly affected via this correlations. However, 
since the relative impact of the first term in Eq. (Q) decreases for increasing a, also the 
effects of the matter density uncertainties should become smaller. 

5 The results for the different measurements 

In this section, we systematically discuss the impact of matter density uncertainties on the 
most interesting measurements at a neutrino factory: sin^ 26*13, the sign of Amg^, and Sqp- 
Although many of these measurements have been individually treated in previous works, 
the goal of this work is to specifically demonstrate where in parameter space and for which 
measurements matter density uncertainties are relevant and important. In addition, we will 
show, where applicable, representative examples of such measurements as functions of the 
amplitude of the matter density uncertainties. We choose the "model of the measured mean 
matter density" with the values Ap* = (no matter density uncertainties), Ap* = 1% ■ p, 
Ap* = 3% ■ p, Ap* = 5% ■ p (the standard value, which we usually propose as a conservative 
choice), and Ap* = 10% ■ p (the most pessimistic choice). 

The sensitivity to sin^ 2^13 

We define the sensitivity to sin^ 26*13 as the largest value of sin^ 26*13, which fits the true 
value sin^ 26*13 = 0. As discussed in the previous section, this choice naturally includes 
the treatment of correlations and degeneracies, which means that any point in the multi- 
dimensional parameter space, which fits the zero rate vector at the chosen confidence level, 
will have a fit value of sin^ 2^13 below the sensitivity limit. Thus, it is guaranteed that 
the experiment will find sin^ 2^13 above the sensitivity limit at the chosen confidence level. 
Any other definition of the sensitivity limit would be more complicated and would involve 
a more sophisticated definition of how to treat degeneracies. For example, one could show 
the sensitive ranges for the different degeneracies separately. Such a definition might be 
sensible for existing experiments, but not to evaluate and compare the potential of future 
experiments in order to optimize them using condensed information. 

As we qualitatively discussed in Sec. El the sensitivity to sin^ 26*13 is limited by the statistics 
in the appearance channels. Therefore, the systematics introduced by the matter density 
uncertainties does not affect it substantially. This can, for example, be seen from comparing 
Fig. 9 with Fig. 10 in Ref. jlU] for NuFact-II. In this impact factor analysis, the matter 
density affects the sin^ 26*i3-precision measurements, but it does not at all affect the sin^ 26*13- 
sensitivity limit. 

The precision of sin^ 2^13 

We define the precision of sin^ 2^13 as the full width error (at the chosen confidence level) 
of the best-fit manifold projected onto the sin^ 2^13-axis. This definition includes correla- 
tions, but it does not include (disconnected) degeneracies, since a sensible definition of the 
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Figure 4: The precision of the measurement o/sin^ 2^13 (full width, best-fit manifold only) 
for NuFact-II and the true value of Sep = as a function of the true value o/sin^26'i3 at 
the 2a confidence level. The different curves correspond to different allowed matter density 
uncertainties Ap* as described in the plot legend, especially the thick curves correspond to 
no matter density uncertainty (light thick curve) and our standard uncertainty Ap* = 5%- p 
(dark thick curve). 

sin^ 2^13-precision including disconnected degeneracies is difficult. In Fig. we show the 
relative error on log(sin^ 2^13) on a double-logarithmic scale. In addition, we only show it 
for one true value of S^p = 0, since otherwise one could not clearly see the effect of the 
matter density uncertainty. However, one should keep in mind that this precision strongly 
depends on the true value of 6cp- A possible representation of this fact is to show the band 
ranging over all true values of the CP phase, such as it is done in Fig. 13 of Ref. |4Uj . 

Figure |3] demonstrates the effects of the matter densities as a function of their amplitude 
Ap* for the LMA-I solution, where the dark thick curve corresponds to our standard choice 
and the light thick curve to no matter density uncertainties, which is often used in other 
studies. Except from the most extreme (and unrealistic) case of Ap* = 10% ■ p, we do not 
find any significant influence of the matter density uncertainties below sin^26'i3 < 10^^. In 
this region, the sin^ 26'13-dependent terms in Eq. (0) are more determined by the statistics 
than the matter density uncertainty. However, at the CHOOZ limit, the matter density 
uncertainties lead to corrections at the percent level of the investigated quantity. It is also 
interesting to note that for Ap* < 1% ■ p the matter density uncertainties are irrelevant. 
Figure 1^ is representative for different values of 6cp than 6cp = 0°. For example, if one 
drew bands using all values of 6cp (such as in Fig. 13 of Ref. [10]), the matter density 
uncertainties would lead to an upward shift of these bands for large values of sin^ 2^13. A 
similar qualitative behavior can also be expected for the LMA-II solution, although the 
general performance with respect to sin^ 2^13 is worse because of the stronger presence of 
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Figure 5: The sensitivity to CP violation for NuFact-II and two different values of 
Amgi as a function of the true values of sin^ 26^13 and 6cp at the 2a confidence level. The 
parameter values within the regions enclosed by the contours refer to the sensitivity to CP 
violation, i.e., the ability to distinguish the respective true value of 6cp given at the vertical 
axis from 6cp G {0,180°} (CP conservation). The different curves correspond to different 
values of Ap* as given in the legend of Fig. ^ especially the thick curves correspond to no 
matter density uncertainty (light thick curve) and our standard uncertainty Ap* = 5% • p 
(dark thick curve). 



correlations and degeneracies. 

The sensitivity to the sign of Amgj^ 

We define that there is a sensitivity to a certain sign of Arrig^, if there is no solution fitting 
the zero rate vector (generated with this specific sign) with the opposite sign of Am|^ at 
the chosen confidence level. This definition implies that the presence of the sgn(Am3^)- 
degeneracy at the chosen confidence level is the most important factor, which destroys 
the sign of the Amg^^-sensitivity. Performing the analysis, it turns out that the matter 
density uncertainty hardly affects the appearance of the degenerate solution. This means 
that matter density uncertainties are not relevant for the sign of the Am|]^-determination, 
at least within the KamLAND-allowed range, even Ap* = 10% ■ p is safe for NuFact-II. 
Given the definition of A at Eq. (^, a compensation of a different sign of Am|^ in A could 
only be achieved by reversing the sign of p, i.e., negative matter densities. However, in 
our analysis, such a reversion would require matter density uncertainties of over 100% ■ p in 
order to allow the unphysical assumption of antimatter in the Earth. 
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The sensitivity to any CP violation 

One can define the sensitivity to any CP violation as the ability of an experiment to dis- 
tinguish any given true CP violating value 5cp ^ {0, 180°} from the CP conserving values 
5cp G {0, 180°}. This definition already implies the treatment of degeneracies: any degener- 
ate solution fitting 5cp ^ {0, 180°} would destroy the sensitivity to CP violation. A special 
case of CP violation covered by this definition is the case of maximal CP violation with the 
true value 5cp ^ {—90°, 90°}, which we will discuss in the next subsection. 

Figure El shows the sensitivity to any CP violation as a function of the true values of sin^ 26iz 
and 5cp for several different values of the matter density uncertainty. The left-hand side 
plot corresponds to the current lower limit of the KamLAND-allowed range, whereas the 
right-hand side plot corresponds to the LMA-I best-fit value. Again, the matter density 
uncertainties are only relevant for sin^ 2^13 > 10~^, as it can be clearly seen in Fig. and 
the measurements are statistics dominated for sin^ 2^13 < 10"'^. Comparing the dark thick 
curves for our standard value Ap* = 5% ■ p between the left- and right-hand side plots at 
the CHOOZ limit indicates that the relative influence of the matter density uncertainties 
decreases with increasing Am^i . In general, the CP performance becomes better for a large 
hierarchy parameter a = Am'^i/Am'^i in Eq. (^, where the CP sensitive terms also become 
large. Thus, the relative impact from the matter density uncertainties in the first term of 
Eq. (0) decreases for larger values of a. In addition, problems with degeneracies are reduced 
by better statistics of the CP sensitive terms. Another important conclusion from Fig. is 
that for Ap* < 1% ■ p the matter density uncertainties are again not relevant. 

The sensitivity to maximal CP violation 

The special case of sensitivity to maximal CP violation 6cp € {—90°, 90°} is represented 
by the horizontal dashed lines in Fig. El This measurement has a much more "binary 
nature" than the sensitivity to a small CP violation, since CP conservation and maximal 
CP violation produce rather distinctive rate vectors for large values of sin^ 26'i3, where matter 
density uncertainties are most important. Since the sensitivity to maximal CP violation is 
close to the optimum of all CP violation sensitivities, it is often used as a representative 
for CP violation discussions. For example, it is shown as a function of the true values of 
sin^ 2^13 and Am\^ in Fig. 17 of Ref. jlHI; a presentation, which would not be possible for 
any CP violation. 

It can be shown that matter density uncertainties are not important for maximal CP vio- 
lation measurements at the 2a confidence level within the (Sa) KamLAND-allowed range 
as long as Ap* < 5% ■ p. This is consistent with Fig. El along the horizontal dashed lines, 
since the left-hand side plot represents the lower limit of the KamLAND-allowed range and 
the relative influence of the matter density uncertainties becomes smaller for larger values 
of Amii. Only for very large matter density uncertainties Ap* = 10% ■ p, they have some 
impact on the very lower end of the KamLAND-allowed range. However, note that this 
statement does not mean that matter density uncertainties in general do not affect the 
(maximal) CP violation measurements. It depends on the used parameter values and it is 
supported by the fact that the KamLAND experiment cut off the lower part of the LMA 
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Figure 6: The precision of the measurement of Scp (full width) for NuFact-II and the 
true value of 6cp = 90° as a function of the true value of sin^ 2^13 at the la confidence 
level. The different curves correspond to different allowed matter density uncertainties Ap* 
as described in the plot legend, especially the thick curves correspond to no matter density 
uncertainty (light thick curve) and our standard uncertainty Ap* = 5%-p (dark thick curve). 



solution. 



The precision of Sqf 

We define the precision of 6cp as the full width error of 5cp on the "CP circle". For a 
given true value of Sep any degenerate solution, which fits the true value, is included in this 
precision. It is sometimes also called the "coverage in Sqp" [101; because it describes how 
much of the CP circle that fits the chosen true value of ^cp- Thus, a precision of 360° (or a 
coverage of 100%) corresponds to no improvement of the knowledge on 6cp. In comparison 
to the CP violation sensitivity, the precision of Sep does not differentiate any special value 
of Sep, such as the CP conserving values and 180°. It describes how much one could learn 
about Scp from a specific experiment. For example, even if an experiment is not sensitive 
to CP violation, because Sqp is too close to CP conservation, it may teach us something 
about Scp and exclude certain regions of the CP circle. In this case, the precision (coverage) 
would be smaller than 360°. 

In Fig. the precision of Sep is shown as a function of the true value of sin^ 2^13 for the 
true value Sep = 90° at the la confidence level. One could also choose a different true value 
of Scp or a different confidence level, but the effects of the matter density uncertainties 
would qualitatively look the same. However, degeneracies are hardly present below the 
la confidence level. Therefore, the numbers on the precision axis should be handled with 
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CP violation discovery reach (3cr) 




1-1 



Figure 7: The 3cr CP violation discovery reach as a function of the true values of sin^ 26'i3 
and 5cp for NuFact-II and JHF-HK from Ref. {jOj - The parameter values within the 
regions enclosed by the contours correspond to the sensitivity to CP violation, i.e., the ability 
to distinguish the respective true value of given at the vertical axis from (5cp € {0, 180°}, 
i.e., CP conservation. The different curves correspond to different values of Ap* as given 
in the figure. For the JHF-HK experiment, matter effects (and matter density uncertainty 
effects) are strongly suppressed by the short baseline L 295 km. For the NuFact-II 
experiment, degeneracies are much more present at the 3a confidence level than at the 2a 
confidence level, as it can be seen by comparison with Fig. [3 (right-hand side plot). 

care, since the 3a errors would be over-proportionally larger not following simple Gaussian 
statistics. As for the case of CP violation, matter density uncertainties become important 
for sin^ 26'i3 > 10~^ and Ap* > 1% ■ p for the same reasons. However, since the CP precision 
measurement includes all possible fit values of 6cp and it is not a binary measurement, the 
effects of the matter density uncertainties can be very large. At the CHOOZ limit, they can 
even affect the precision by a factor of two or more. 

6 Summary and conclusions 

Matter density uncertainties up to 5% can be present on the PREM matter density profile. 
It is well known that they can influence many of the measurements at a future neutrino 
factory, such as the CP precision measurements. In this work, we have first of all discussed 
different approaches to model the matter density uncertainties and how they are related to 
each other. We have concluded that using the average matter density with an uncertainty 
of about Ap* = 3% - 5% is an appropriate conservative choice for a complete statistical 
neutrino factory simulation. In this approach, the average matter density is measured by 
the neutrino factory as an independent parameter within the externally given precision Ap*. 



16 



It has several advantages: 

• It is fast enough for a complete statistical simulation including systematics, correla- 
tions, and degeneracies. 

• It allows correlations between the matter density and the neutrino oscillation param- 
eters. 

• It can use the external information from other statistical models or a specific baseline 
to improve the knowledge on the input parameter Ap*. 

• It can even simulate the matter density profile effect for not too long baselines (L < 
5 000 km). 

With this model, we have systematically discussed how the matter density uncertainties 
affect the most important measurements at a large neutrino factory at a baseline of L = 
3 000 km. We have found that some, but not all, measurements suffer from the matter 
density uncertainties. In addition, they are only relevant in certain regions of the parameter 
space. In detail, we have found that 

For the sin^ 20i3-sensitivity matter density uncertainties are not important, because 
this measurement is constrained by the statistics in the appearance channels and 
not by the systematics coming from the matter density uncertainties. 

For the sin^ 20i3-precision matter density uncertainties are only relevant for sin^ 2^13 > 
10~^ because of the statistics domination below that region. 

For the sign of Arrig-^- sensitivity matter density uncertainties are irrelevant because of 
the very distinctive solutions of this binary measurement, which could only be mixed 
up by the matter density uncertainties by using negative matter densities. 

For the sensitivity to any CP violation matter density uncertainties again become im- 
portant for sin^ 2^13 > 10~^. However, their relative influence is reduced for large 
values of Aml^. 

For the sensitivity to maximal CP violation matter density uncertainties are not sig- 
nificantly important within the KamLAND-allowed range at the 2a confidence level. 

For the precision of <5cp matter density uncertainties are not important for sin^ 2^13 < 
10~^. However, for sin^ 29i3 > 10~^ they make the precision worse by a factor of two 
(close to the CHOOZ hmit). 

We conclude that three conditions are necessary for the matter density uncertainties to be- 
come relevant at the considered neutrino factory analysis: First, only some of the measure- 
ments are influenced by matter density uncertainties, such as the sin^ 2^13- and 5cp-precision 
measurements. Second, matter density uncertainties are only important for sin^ 2^^i3 > 10^^. 
Third, matter density uncertainties only sensibly affect the measurements for uncertainty 
amplitudes Ap* > 1% • p, as it can be observed in any of the figures so far. 
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Since superbeams or superbeam upgrades are basically sensitive down to sin^ 26*13 ~ 10^^ 
for the discussed measurements, competitiveness of neutrino factories with superbeams is an 
interesting issue. Superbeams are often proposed with shorter baselines and they are using 
different (non-resonant) neutrino energies, which means that they are much less affected by 
matter density uncertainties. This is one of the reasons why they are often much better than 
the neutrino factories for sin^2^i3 > 10~^. For illustration, the 3a CP violation discovery 
potential is shown in Fig. [7| for NuFact-II (with Ap* = 5% ■ p and Ap* < 1% ■ p) and the 
JHF to Hyper-Kamiokande superbeam upgrade JHF-HK from Ref. jlU] (as it is proposed 
in Ref. jTTj, using two years of neutrino running and six years of antineutrino running). 
Since this superbeam upgrade is proposed with a quite short baseline L ~ 295 km, it is 
hardly affected by matter effects and matter density uncertainty effects. The figure clearly 
demonstrates that especially for sin^ 2^13 > 10~^ a knowledge on the matter density profile 
with about 1% precision is necessary for the competitiveness of neutrino factories with 
superbeams. Therefore, it would be interesting to know with what costs such a precision 
could be achieved in geophysics along a specific baseline. However, if the superbeams do not 
find sin^ 2^13 > 0, then the neutrino factories will for sin^ 2^13 < 10~^ be almost unaffected 
by matter density uncertainties in the post-superbeam era. 
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